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Peak detection is one of the most important time-domain functions performed in signal monitoring. Peak detection is the 
process of finding the locations and amplitudes of local maxima and minima In a signal that satisfies certain properties. 
These properties can be simple or complex. For example, requiring that a peak exceeds a certain threshold value is a 
simple property. However, requiring that a peak s shape resembles that of a prototype peak is a complex property. 

Peak detection is important in many applications, such as chemistry, biology, and music. Scientists and engineers who 
use analysis techniques such as spectroscopy, chromatography, and tone detection often use peak detection methods 
specific to those analysis techniques. However, this document describes a general method that applies to a variety of 
signal types. This is the method used in LabVIEW and Measurement Studio for peak detection functions. 



In some applications, you do not need to know the exact peak amplitudes and locations, rather you need to know the 
number or general locations of peaks. In this case, use a threshold peak detection function, such as the Threshold Peak 
Detector VI in LabVIEW. 

Figure 1 shows the Threshold Peak Detector Viand the Vl's inputs and outputs. The VI scans the Input sequence X, 
searches for valid peaks, and keeps track of the indices of the beginnings of the peaks and the total number of peaks 
found. A peak is considered valid if it has the following characteristics: 

The elements of X begin below threshold, exceed threshold at some index, and then return to a value below 

threshold 

The number of successive elements that exceed threshold is greater than or equal to width. 




Figure 1 . Threshold Peak Detector VI 



This Vi does not identify the locations or the amplitudes of peaks with great accuracy, but the VI does give an idea of 
where and how often a signal crosses above a certain threshold value. 




Some applications require more robust and accurate pBak detection algorithms. The rest of this document focuses on 
uses of advanced peak detection functions, tips to keep in mind and pitfalls to avoid while using them, and methods for 
ensuring that your peak detection measurements are accurate and useful. 

The following sections focus mainly on peaks. However, except where noted, the same information can be used for 
finding valleys or local minima. 



Figure 3 shows the Peak Detector VI and the Vl's inputs and outputs. Figure 4 shows the equivalent function prototype in 
the Advanced Analysis library of LabWindows/CVI; ComponentWorks and ComponentWorks++ contain similar 
interfaces for this function. 




Figure 3. Peak Detector' 




Figure 4. PeakDetector function prototype for LabWindows/CVI 



Notice that it takes handles to the locations, amplitudes, and second derivatives arrays. For descriptions of each of the 
inputs and outputs for the Peak Detector VI orthe PeakDetector function, refer to the LabVIEW or Measurement Studio 



http://zone.ni.corn/dewone/cda/tut/p/id/3770»tocl 
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LabVIEW also contains a Waveform Peak Detection VI as shown in Figure 5. 




Figure 5. Waveform Peak Detection Vi 



The Waveform Peak Detection VI operates like the array-based Peak Detector VI. The difference is that this Vl's input is a 
waveform data type, and the VI has error cluster Input and output terminals. Locations displays the output array of the 
peaks or valleys, which is still in terms of the indices of the input waveform. For example, if one element of Locations is 
100, that means that there is a peak or valley located at index 1 00 in the data array of the input waveform. Figure 6 

vs you a method for determining the times at which peaks or valleys occur. The following equation locates the peaks 



and valleys: 



Time Locatlonsp] = tO + dt*Locatlons[i] 




Initialize (T)llC 



Figure 6. Using the Waveform Peak Detection VI to i 



the times at which peaks or valleys oc 



usefulness of peak detection 

1 . The function can process many blocks of data that are part of the same signal. By correctly using the Initialize (T) and 
end of data (T) inputs, you can use the peak detection function to analyze a signal that has been broken up into several 
data blocks. You can also acquire a continuous signal and process pieces of the signal as they become available. The 
VI finds the peak locations in each block, relative to the previously analyzed blocks. For example, to process a signal 
acquired in five consecutive blocks, you can use the following pseudocode algorithm: 

for i - 1 to 5 



Initialize = False; 

if (i == 5) 
EndOfData - True 
else 

EndOfData = False; 

Set polarity (peaks or valleys), width, threshold 
Call PeakDectector function 

Copy the output values to different variables so they will n 



The same algorithm in LabVIEW might look like the VI in Figure 7. 




lg the Waveform Peak Detection VI to process a signal that is broken into several blocks 



This multiple-block feature al 



-ailable. The data bl 



The function retains history information; it uses a history buffer to retain a certain number of data points from the previous 
data block. This feature allows the function to correctly locate peaks and valleys that are close to the boundary between 
blocks. However, you must set end of data (T) and Initialize (T) when you have finish analyzing one signal and are 
starting a new one. Otherwise the function will view the two blocks as part of the same continuous signal and will 
incorrectly locate peaks or valleys near the end of the last block of the previous signal, or near the beginning of the first 
block of the new signal. 

3. The peak location function gives peak locations at fractional indices. It uses the quadratic fitalgorithm and returns the 
peak locations as floating point numbers, not as integer index values Therefore, the peak locations and amplitudes 
usually do not correspond to actual data points in the sampled input signal. 

This feature is an advantage of the algorithm because it effectively interpolates between the data points while finding 
peaks and valleys. The function can therefore measure peaks that have a greater amplitude than any data points near 
the peak. This interpolation provides a good indication of the true value of the peak in the original analog signal. 

4. The function allows implicit noise reduction while finding the peaks. Using the width parameter in some cases can 
effectively reduce the noise in the input signal when finding the peaks. The minimum value is three; using this value 
results in no noise reduction. Using a width value larger than three implicitly smooths the data. This feature is useful in 
some applications. However, you must ensure that you use large width only on noisy data. Y 
the peak locations and amplitude results are reasonable. 



whether a peak or valley is present. 

If width = 3, then the fit is exact, meaning the parabola will actually pass through each of the three points. If width is 
greater than three, then a least-squares fit is performed. This process will smooth high-frequency noise if the width is 
sufficiently large. 

For each set of points, the algorithm performs the least-squares quadratic fit, and then performs a series of tests on the 
coefficients to see whether they meet the criteria for a peak. The function checks whether each parabola is at a local 
maximum, determines the sign of the quadratic coefficient, which Indicates the parabola s concavity, and finally checks 
that the peak is above the designated threshold. 

• 



To use the peak detection function correctly, you need to have some prior knowledge about your signal. The following 
are some important issues to consider when specifying the input parameters: 
Is the data a time-domain or a frequency-domain signal? 

If the data is a frequency domain signal and contains one or more well-defined frequency components, use the Power 
& Frequency Estimate VI in LabVIEW or the PowerFrequencyEsbmate function in Measurement Studio. These 
functions allow you to get more accurate information aboutthe exact frequencies of the peaks and their 
corresponding energies. LabVIEW 6.0 also includes the Extract Single Tone Information VI, which uses a curve-fitting 
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method to precisely identify the amplitude and frequency of the largest frequency component of a signal. 
Are all the peaks that you are looking for at roughly the same amplitude? 

If so, then set one threshold value for analyzing all the data. However, if you expect to have peaks at many different 
amplitudes, break up the data and use a different threshold value for separate data blocks. 

Is your data periodic? 

If so, then searching for more than a few periods of the data for peaks is an inefficient use of processing time. Often 
there is some noise or other distortion present in the data. In this case, average many periods of the original time 
signal to get one or a few averaged periods, then pass these to the peak detection function. 
Is the data adequately sampled? 

Inadequate sampling of data can result In Inaccurate values for the locations and amplitudes of peaks, and the non- 
detection of valid peaks. Despite the fact that peak detection is essentially a time-domain operation, the sampled 
signal must still satisfy the Nyquist sampling theorem the sampling rate must be at least twice the largest frequency 
component In the signal. However, useful digital representation of a signal typically requires a sampling rate between 
five and ten times the largest frequency component The front end of the data acquisition system should contain an 
analog anti-aliasing filter that removes or strongly attenuates components above Nyquist sampling rate, or half the 
sampling rate. 
Is your data noisy? 

This is a common problem that you must deal with very carefully in peak detection. Clearly, high-frequency noise 
results in the detection of a large number of peaks, but typically only a few of these will actually be of interest In these 
cases, increase the width parameter to implicitly smooth the data for finding the peaks. Notice that since this process 
tends to remove high-frequency spikes from the data, Increasing width tends to decrease the amplitudes of peaks and 
increase the amplitudes of valleys. You can also use explicit smoothing or interpolation techniques to effectively 
reduce noise in your data. 



Smoothing can cause problems with the measurements if used incorrectly. You need to conduct some initial testing to 
determine the optimum width. This means you need to determine the number of data points that will remove a sufficient 
amount of noise without removing significant features of the original signal. However, since you cannot explicitly see the 
smoothed data, it is difficult to determine the optimum width through visual verification. 

In most cases, it is preferable to smooth or process the data before applying the peak detection function. In that case, 
use a width of three, which tells the peak detection function to process exactly the signal you passed to it, without any 
smoothing. This method gives you more control over the smoothing and processing of the original signal. The peak 
detection function then processes only the signal that you pass to it 



Figure 8b shows the same signal analyzed by the peak detection function using width = 29, a relatively large value 
function detects only three peaks in this case. These peaks may be correct, but it is difficult to determine whether thi 
peak information is really useful. 




Figure 8a. A noisy waveform after being passed through the peak d< 

;kdots mark the detected peaks; most of the marked peaks are not really of interest, butar 
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Figure 8b. The same signal as in Fig. 8a using width = 29 

turned by the function. It is difficult to tell whether the dots represent 
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There are several common methods of interpolating discrete-time signals, A simple method is linear interpolation. If you 
have an interpolating factor of q, then q-1 points are inserted between each original data point. With linear interpolation, 
the added points lie on the line through the original data points on either side. While these results are not an accurate 
representation of the original analog signal, they can help the peak detection algorithm to detect all of the valid peaks. 

Another method of interpolation is to interleave q-1 zeros between each of the original data points, then execute a 
lowpass filter. The theory behind this method is beyond the scope of this document; however, many texts on discrete- 
time signal processing contain an explanation of this method. For more information on this method, refer to the 
References section at the end of this document. 

A certain trade-off comes with interpolation. Interpolating the digitized signal tends to place the found peaks closer to 
actual points in the interpolated signal. However, depending on the nature of the anBlog signal, these may or may not be 
closer to the real peaks than those detected using the uninterpolated digital signal. So the trade-off is between finding all 
valid peaks and getting more accurate data for the peak locations and amplitudes. 



You now know how to use the peak detection functions in LabVIEW and Measurement Studio to find the locations and 
amplitudes of peaks and valleys in your signals. The powerful features of these functions lend themselves to convenient 
and accurate signal analysis. However, as described above, it is important to understand the significance of the input 
parameters when using the functions. Furthermore, to use these functions effectively, you need to have an 
understanding of the nature of the input signal before using the peak detector functions. 




The resources linked below contain more information about the theory behind methods of digital signal processing, 
specifically digital filtering, and interpolation. For valuable information about frequency-domain analysis, digital filtering, 
and using LabVIEWand Measurement Studio to apply the ideas presented in this document, refer to the resources 
linked below. 

You can order the following books through National Instruments Books and Publications linked below: 

Chugani, Mahesh L, Samant Abhay R., and Cerna, Michael. LabVlBW Signal Processing, Prentice Hall, 1998. 
Haykin, Simon and Van Veen, Barry, Signals and Systems. John Wiley & Sons, Inc.. 1998. 

V. and Schafer, Ronald W. Discrete Time Signal Processing. Prentice Hall, 1999. 



Very Helpful 

The above description of Peak Detection gives a good basic explanation of the way the function detects peaks and 
how the parametres affect the outcomes. 

-jaspal.boparai@gaiam.co.uk - Jan 23, 2009 



no explanation about least squares approximation 

in the article there's no explanation about the least squqres approximation method, i believe, u should include a short 
description of it. 

- Nathan Minaev, Tel Aviv University, Israel, mina8van@posttau.ac.il - Oct 27, 2005 



Legal 

This tutorial (this tutorial-) was developed by National Instruments ("Nl"). Although technical support of this tutorial 
may be made available by National Instruments, the content in this tutorial may not be completely tested and venfied 
and Nl does not guarantee its quality in any way or that Nl will continue to support this content with each new revision 
of related products and drivers. THIS TUTORIAL IS PROVIDED "AS IS" WITHOUT WARRANTY OF ANY KIND AND 
SUBJECT TO CERTAIN RESTRICTIONS AS MORE SPECIFICALLY SET FORTH IN NI.COM'S TERMS OF USE 
(http://ni.co m/legal/termsofuse.'unitedstates/us/). 
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